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CRITICAL SCALING FOR THE SIMPLE SIS STOCHASTIC 

EPIDEMIC 

R. G. DOLGOARSHINNYKH AND STEVEN P. LALLEY 


Abstract. We exhibit a scaling law for the critical SIS stochastic epi¬ 
demic: If at time 0 the population consists of y/N infected and N — s/~N 
susceptible individuals, then when time and number currently infected 
are both scaled by \fN, the resulting process converges, as N —> oo, to 
a diffusion process related to the Feller diffusion by a change of drift. 
As a consequence, the rescaled size of the epidemic has a limit law that 
coincides with that of a first-passage time for the standard Ornstein- 
Uhlenbeck process. These results are the analogues for the SIS epidemic 
of results of Martin-Lof |§] for the simple SIR epidemic. 


1. Introduction 

Among the most thoroughly studied stochastic epidemic models are the 
simple SIR, and SIS epidemics (see Weiss and Dishon m for the origin of 
the SIS model); and among the many problems associated with these models 
perhaps the most basic and most interesting have to do with the duration 
and size of the epidemic. When the epidemic is either sub- or supercritical, 
the large-population behavior of the duration and size is reasonably well- 
understood: in the subcritical case, the epidemic is stochastically dominated 
by a subcritical Galton-Watson process (see below), and in the supercritical 
case the epidemic may become endemic (see Norden m and Kryscio and 
Lefevre 0). For critical epidemics, large-population asymptotics are more 
delicate. In 1998, Martin-Lof |5] (see also Aldous li) showed that the size 
S = Sn has an interesting and nontrivial asymptotic behavior as the popu¬ 
lation size N —> oo: If the number Ao of individuals initially infected is of 
order 6A * 1 / 3 , then the size <Sjv (defined to be the sum of the total infection 
times over the whole population) has a limit distribution 

(1) S n /N 2 / 3 

where T fe * is the first passage time of W(t ) + t 2 /2 to the level b, and W(t) 
is a standard Wiener process. Furthermore, 1/3 is the critical exponent, 
in that the quadratic drift is not felt if Ao is much smaller than A 1 / 3 : If 
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X(0) rs_/ bN a for some a < 1/3 then S]\r/N 2a converges in law to the first 
passage time of W(t) to the level b. (Martin-Lof neither proved nor stated 
this, but it can be deduced from his methods.) 

The purpose of this note is to establish an analogous scaling law for the 
critical SIS epidemic: If the number Xq of individuals initially infected is of 
order b\/N, then as N —> oo, 

(2) S n /N Tb , 

where T b is the time of first passage to 0 by a standard Ornstein-Uhlenbeck 
process started at b. Our approach has a rather different character than 
those of Aldous and Martin-Lof. We shall establish that the SIS epidemic 
process itself, suitably rescaled, converges in law to a diffusion process we 
dub the attenuated Feller diffusion. This has a law absolutely continuous 
with respect to, but not equal to, that of the standard Feller diffusion. (See 
equation © below for the stochastic differential equation governing the 
Feller diffusion, and © for the attenuated Feller diffusion.) Furthermore, 
we will show that 1/2 is the critical exponent, in the following sense: If 
X(0) rs_/ bN a for some a < 1/2 then the rescaled SIS process converges 
in law to a standard Feller diffusion with no drift. It will follow that the 
duration of the epidemic, rescaled by yW(0), converges in law as N —> oo 
to the first passage time to 0 of the corresponding Feller diffusion. 

Our analysis will show that the critical scaling window for the transmis¬ 
sivity parameter /3 (see below) is on the order 1 /y/N. This scaling window 
has also been observed — but in a different context — by Nasell m, nn 
Nasell shows that the quasistationary distribution of the SIS epidemic un¬ 
dergoes a scaling transition when the transmissivity parameter varies from 
below 1 — 0(N _1 / 2 ) to above 1 + 0(N ~ 1//2 ). This phenomenon does not 
seem to be directly linked to the critical scaling in our Theorem El 

2. The SIS Epidemic and its Branching Envelope 

2.1. SIS Model. The SIS epidemic is a continuous-time birth-death Markov 
chain Xf = X(t) on the state space [N] := {0,1, 2,... , N} whose infinitesi¬ 
mal transition probabilities are as follows: 

(3) P{X(t + St) = x + 11 X(t) = x} = f3x( 1 — x/N)5t + o(xSt), 

P{X(t + St) = x — 11 X(t) = x} = xdt + o(x6t ), 

P{X(t + St) = x | X(t) = x} = 1 — f3x( 1 — x/N)5t — xdt + o(xSt)). 

These describe a population of N individuals in which Xf are infected and 
the remainder N — Xf are susceptible to infection at time t. Infected indi¬ 
viduals recover at rate 1, after which they once again become susceptible to 
(re-)infection, and susceptible individuals become infected at rate (3X t /N 
proportional to the number of infected individuals in the population. The 
epidemic ends at the first time T = Tjy = t when Xf = 0 (note that state 
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0 is absorbing). The epidemic is said to be critical when 0 = 1, and nearly 
critical when 0=1 + A /\/N. 

2.2. Branching Envelope. When the number of individuals infected is 
small compared to the population size, the epidemic evolves approximately 
as a continuous-time branching process Z(t ) = Zt with infinitesimal transi¬ 
tion probabilities 

(4) P{Z{t + 5t) = z + 1 1 Z(t) = z} = 0z5t + o(z5t), 

P{Z(t + 5t) = z — 11 Z(t) = z} = zSt + o(z8t), 

P{Z(t + 8t) = z | Z(t) = z} = 1 — (0 + 1 )z5t + o(zSt). 

We shall refer to this process as the branching envelope of the SIS process. 
Observe that the death rate x is the same as for the SIS epidemic, but the 
the birth rate 0x dominates the birth rate 0x(l — x/N) of the SIS process; 
the difference 0x 2 /N will be called the attenuation or attrition rate. It is 
possible, by a standard construction, to build the SIS process X{t) and its 
branching envelope Z(t ) on the same probability space in such a way that 
X (0) = Z( 0) and X(t) < Z(t ) for all t > 0. Thus, the size and duration 
of the SIS epidemic are stochastically dominated by the total progeny and 
extinction time of the branching envelope. 

2.3. Critical Scaling for the Branching Envelope. It was proved by 
Feller [I] that a critical branching process, when properly renormalized, 
behaves approximately as a Feller diffusion with drift AY), that is, a solution 
of the stochastic differential equation 

(5) dY t = AY) dt + y/YtdW t 

where Wt is a standard Wiener process. (Equivalently, the Feller diffusion 
with drift parameter A may be described as the diffusion process on [0, oo) 
with infinitesimal generator Q x = \xd x + xd 2 x /2.) Feller’s theorem (see 
Jirina {Bj and Lindvall 8] for the proof) asserts that if Z m (t ) is a sequence 
of branching processes satisfying Q with 0 = 0 m = 1 + X/m and with 
Z rn { 0) ~ bm for some b > 0 and A 6 R then 

(6) Z m (mt)/m —^ Y) 

where Y) is the Feller diffusion with drift parameter A and initial value 
Y 0 = b. 

2.4. Critical Scaling for the SIS Process. Because the branching enve¬ 
lope stochastically dominates the SIS process, the scaling law © limits the 
duration and growth of the critical and near-critical SIS process. Since time 
is scaled by the factor m, where Z m ( 0) ~ bm, it follows that the correspond¬ 
ing SIS started with X (0) ~ bN a infected individuals cannot have duration 
longer than Op(N a ) time units. Consequently, we should expect that if the 
attenuation rate, divided by the scale factor N a and integrated to time N a , 
is op(l) then the limiting behavior of the rescaled SIS process X(N a t)/N a 
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should be no different from that of the branching envelope Z(N a t)/N a . An 
easy calculation shows that this will be the case when a < 1/2. When 
a = 1/2, the accumulated attrition over the duration of the branching enve¬ 
lope will be on the same order of magnitude as the fluctuations, and so the 
rescaled SIS process should have a genuinely different asymptotic behavior 
from the branching envelope. Our main result makes this precise: 

Theorem 1. Assume that process X(t) = X N (t) has infinitesimal transition 
probabilities W- If for some constants a < 1/2 and b > 0 the number of 
individuals initially infected satisfies X N (0) ~ bN a , and if the birth rate n 
satisfies /3 = fix = 1 + X/N a , then as N —> oo, 

(7) X N (N a t)/N a Y t 
where 

(a) if a < 1/2 then Yt is a Feller diffusion with drift X and initial state 
Y 0 = b; 

(b) if a = 1/2 then Y t is an “attenuated” Feller diffusion with drift A 
and initial state Yq = b, that is, Yt is a solution to the stochastic 
differential equation 

(8) dY t = (A Y t - Y t 2 ) dt + y/YtdWt 
with Wt a standard Wiener process. 

Note that the “attenuation” term —Y 2 in the drift of the limiting process 
© can be guessed from the form of the attrition fix 2 /N. The proof of 
Theorem HI will be given in section 0] below. 


3. Size of the Epidemic 


The size of an epidemic can be defined as the total number £ of new 
infections during its entire course. Alternatively, it can be defined as the 
total infection time summed over the whole population: 

rT 


(9) 


S 


Sn 


X t dt. 


Jo 

Although the two definitions are not the same, it can be shown that the 
two quantities have the same asymptotic behavior for large N, that is, £ at ~ 
Sn- (This follows from the fact that the length of the infection periods for 
infected individuals are i.i.d. unit exponential r.v.s.) Because the integral © 
is a continuous functional of the path X t (relative to the Skorohod topology), 
Theorem |T| implies that if X(0) ~ b\fN and fi = 1 + X/y/N then 


v r (0) 

(10) S n /N-X> / Yt dt 

Jo 

where Yt is the attenuated Feller diffusion Q with initial state Yq = b and 
tq is the first passage time to 0 by Y t . 
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By an odd bit of luck, the instantaneous rate Y) dt at which infection time 
accrues coincides with the rate of change in accumulated quadratic variation 
of the semimartingale Yt. (Note: In fact this is really no accident, but rather 
an artifact of the fundamental connection between Galton-Watson processes 
and random walks via the “depth-first search” algorithm. See jl] for more 
on this.) This suggests making the natural time change to the diffusion Yt 
so as to make the instantaneous quadratic variation constant. The new time 
scale s = s(t) and the old t are related by 

(11) ds = Y t dt, 

and so f Y t dt = f ds is the limit of the rescaled epidemic sizes Sjsr/N. The 
time-changed process V s = Y t ^ satisfies the stochastic differential equation 

(12) dV s = (A - V s ) ds + dW s 

where W s is again a standard Wiener process. Setting U s = V s — A, one ob¬ 
tains the stochastic differential equation for the standard Ornstein-Uhlenbeck 
process: 

(13) dU a = -U s ds + dW s 
This proves the following. 

Corollary 1. If X (0) ~ b\^N and (3 = 1 + A /yfN then 

(14) Sn /N^t( 6-A;-A) 

where r(x; y) is the time of first passage to y by a standard Ornstein- 
Uhlenbeck process started at x. 

The Laplace transforms of the distributions of t{x\ y) can be expressed 
in terms of parabolic cylinder (Weber) functions: see Darling and Siegert 
J2j. These do not invert easily. However, in the special case A = 0 (the case 
corresponding to the critical SIS epidemic!), the distribution of r(6; 0) has 
a simple closed form: 

(15) P{r(b ; 0) > s} = P b {U s > 0} - P b {U s < 0}. 

This can be obtained from a reflection principle, using the symmetry of the 
Ornstein-Uhlenbeck process about the origin. 

4. Proof of Theorem □ 

We prove Theorem Q using the weak machinery developed in Ethier and 
Kurtz |31, which reduces the problem to checking convergence, in an appro¬ 
priate sense, of infinitesimal generators. Write 

(16) 


Y t N = X N {N a t)/N a 
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for the rescaled epidemic process, and denote by Ey the corresponding 
expectation operator under the initial condition Y^ 1 ' = y. For For / E 
C[ 0, oo), define 


(17) 

(18) 
(19) 


and 

h—> 0 h 

Gf(y) = (ty)-^f(y) + y--^(y) tfa<i/2, 
Gf(y) = (Ay - y 2 ) ■ jf(y) + y ■ ^-£(y) if a = 1/2 


dy 


dy 2 


By Corollary 1.2, sec. 11.2 of 2j, the operator Q restricted to C[0,oo) PI 
(7 2 (0, oo) generates a Feller semigroup on C7[0, oo), and Proposition 3.3, 
sec. 1.5 implies that C'^°[0,oo) is a core for the generator. (An easy calcu¬ 
lation, which we omit, shows that 0 is an exit boundary and oo is a natural 
boundary in both cases.) Moreover, the Markov processes determined by 
these Feller semigroups can be constructed so as to satisfy the stochastic 
differential equations © and ®, respectively. By Theorem 2.5 on page 167 
and Theorem 6.1 on page 28, to prove convergence 0 it is enough to show 
that for each / in the core of Q the generators converge in the sense of the 
following lemma. 


Lemma 1. Let f E C' ( ?°[(l, oo) then 

(20) lim sup | Q N f{y) - Gf(y\ = 0. 

°V[iV]/V 

Proof. Consider first the case a = 1/2. The first step is to calculate Q N f 
for / E oo). Using the infinitesimal transition probabilities (121 1) , we 

have (with x = yN 1//2 and h = tN 1(/2 ) 


= /(1 + A 1 / 2 ) - f(y)] X f(l + AA~ 1 / 2 )yA 1 / 2 (l - yN^/^hN 1 / 2 


+ 


/(1 - A" 1 / 2 ) - f{y) 1 x [yA 1 / 2 • hN 1 / 2 ] + o{Nhy) 


/(l + N- 1 ' 2 ) - /(y)] x [(yN 1 ' 2 + Ay - y 2 - Ay 2 A" 1 / 2 )hA 1 / 2 

+ [/(1 - A" 1 / 2 ) - /(y)] x [yN 1 / 2 ■ hVN] + o(Nhy). 

The error term o(Nhy) is uniform in y, because / is assumed to have com¬ 
pact support. Taking the limit of this expression as h —*• 0 yields 

G N f(y) = (/(1 + iV- 1/2 ) - f(y)) X A 1 / 2 (Ay - y 2 - Ay 2 A" 1 / 2 ) 


+ 


(/(l _ A- 1 / 2 ) - /(y)) - (/(y) - /(y - A -1 / 2 )) 


x Ay. 
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Since / E C^°[0, oo), there exists a constant C > 0 so that / and all its 
partial derivatives vanish for y > C . Therefore, uniformly in all y E \N]/yN, 

v™J N f{y) = x (Ay - y 2 ) + |^(y) x y = Gf(y). 

A similar calculation establishes convergence of generators when a < 

1 / 2 . □ 


5. The SIR Epidemic Revisited 


The continuous-time SIR epidemic differs from the SIS epidemic in that 
individuals may only be infected once: upon recovery, individuals are ef¬ 
fectively removed from the population. Thus, the state at any time t is 
determined by two variables, the number currently infected (I(t) = I N (t )) 
and the number removed ( R(t ) = R N (t)). These take valued in the set of 
nonnegative integer pairs (i, r ) such that 0 < % + r < N, where N is the 
(original) population size. The instantaneous transition rates are as follows: 

(21) (i, r) *—> (i — 1, r + 1) at rate idt ; 

(i,r) I—>■ (i + l,r) at rate f3i(N — i — r)dt/N. 

All states (i, r ) with i = 0 are absorbing: the epidemic ends at the first time 
one of these states is visited. 

As for the SIS epidemic, if the numbers of infected and removed individ¬ 
uals are small compared to the total population size N, then the second 
transition rate in m reduces to (3idt , and so the process I(t) evolves ap- 
proxiamtely as the branching process (J3}. Therefore, by the same logic as in 
sec. I2~il the limiting behavior of the epidemic can be deduced by examina¬ 
tion of the accumulated attrition over the duration of the branching process. 
The result is as follows. 


Theorem 2. Assume that {I N (t), R N (t)) has instantaneous transition rates 
m, and assume that R N ( 0) = 0. If for some a < 1/3 and b > 0 the number 
I N {0) 0 j individuals initially infected satisfies I N ( 0) ~ bN a , and if the birth 
rate (3 satisfies (3=1 + \/N a , then as N —> oo, 


( 22 ) 



(m 


where (i) if a < 1/3, the limit process ( I(t),R(t )) satisfies 


(23) dl(t) = A I{t) dt + sjJ{i ) dW t 

dR(t ) = I(t) dt ; 


and (ii) if a = 1/3, 

(24) dl(t) = (XI(t) - I(t)R(t)) dt + yfUfy dW t 

dR(t) = I(t ) dt. 
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Martin-Lof’s result HJ can be easily recovered from Theorem [21 by the 
same device as used in sec.EJabove. Define the new time scale s by ds = It dt 
and corresponding time-changed process dJ(s) = dl(t). Then the total size 
of the epidemic is just the integral f ds up to the time of first passage to 0 
by J(s). But J(s) is just the Wiener process with a quadratic drift, and so 
IH) follows. 

Theorem [21 can be proved either by martingale methods or by use of the 
Ethier-Kurtz machinery. The latter approach is mildly complicated by the 
fact that the generator 

(25) Q = (A i — ir)di + Vida + id r 

is not elliptic, but rather parabolic, and singular along the i = 0 axis. 
The singularity at i = 0 can be handled by truncating the state space: To 
prove (123), it suffices to prove weak convergence for the processes (I N (t A 
r £ ),R N (t A r e )), where t £ is the time of first passage to the level i = e. 
Nonellipticity of the generator may be handled by using standard existence 
results from the theory of parabolic PDE (see (Sj) to verify the hypotheses of 
the Hille-Yosida theorem (Th. 2.2 of [2]). Weak convergence of the truncated 
processes may then be proved by checking convergence of generators; this is 
another routine calculation similar to that carried out for the SIS epidemic 
in sec. 0 above. 
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